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ABSTRACT 

We use photometric and spectroscopic infrared observations obtained with the 
Spitzer Space Telescope of 12 radio-loud active galactic nuclei (AGN) to investigate 
the dust geometry. Our approach is to look at the change of the infrared spectral 
energy distribution (SED) and the strength of the 10 /im silicate feature with jet 
viewing angle. We find that (i) a combination of three or four blackbodies fits well the 
infrared SED; (ii) the sources viewed closer to the jet axis appear to have stronger 
warm {^ 300 — 800 K) and cold (~ 150 — 250 K) dust emissions relative to the 
hot component; and (iii) the silicate features are always in emission and strongly 
redshifted. We test clumpy torus models and find that (i) they approximate well the 
mid- infrared part of the SED, but significantly underpredict the fluxes at both near- 
and far-infrared wavelengths; (ii) they can constrain the dust composition (in our 
case to that of the standard interstellar medium); (iii) they require relatively large 
(~ 10% — 20% the speed of light) redward displacements; and (iv) they give robust 
total mass estimates, but are insensitive to the assumed geometry. 
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1 INTRODUCTION 

Some active galactic nuclei (AGN) exhibit in their optical 
spectra both broad (~ 1% — 5% the speed of light) and 
narrow emission lines (type 1 AGN) and some only nar- 
row (type 2 AGN). The existence of apparently two distinct 
classes of AGN has been explained within unified schemes 
by orientation effects: an optically thick, dusty torus (or 
warped disc) located outside the accretion disc obscures the 
broad e mission line region for some lines of sight (see re- 
views bv lLawrencelll987l : lAntonuccilll993l : IUrrv fc Padovanil 
1 19951 ). The strongest observational evidence interpreted in 
favour of such a torus are broad emission lines in the 
polarized, scattered liRh t of numerous typ e 2 AGN (e.g.. 



Antonucci fc Milleij[l985l: ICohen et all 19991 : iLumsden etal] 



200 ll ) and significant infrared emission in the continuua 
of most AGN. However, details of the physical state of 
this torus, its precise geometry and location remain uncon- 
strained to this day. 

The putative torus will absorb a significant fraction of 
the nuclear (accretion) luminosity and reradiate it strongly 
in the infrared. The dust radiative trans fer problem was first 
invest igated for a toroidal g eomet r y by Krolik fc BegelmanI 
l|l988h and iPier fc Kroli"^ (|l992l . Il993l '). Because of the 
difficulties in modelling a clumpy medium they assumed a 
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uniform density distribution, although they noted that the 
dust must be concentrated in clouds to protect the grains. 
Their models reproduced the gross features of the observed 
spectral energy distributions (SEDs), but some major 
problems remained: (i) the torus was predicted to emit 
anisotropically; (ii) the generated emission at far-infrared 
wavelengths was insufficient; (iii) early observations of 
type 1 AGN did not sho w the predicted s t rong 10 fj,m 
silica t e emission featur es (|Roche et al.l Il99ll : IClavel et al.l 
l200d : ISpoon et al.ll2002r ): and (iv ) the required geometrical 
thick ness was difficult to support (|Krolikll2007l : IShi fc KrolikI 
l200a ). These problems persisted in all the following studies, 
which also employe d cont in uous density distributions (e.g., 
Granato & Dancso 1994: Efstathio u fc Rowan- RobinsonI 
1995; van Bcmmcl & DuUcmond 2003). 

More recently, iNenkova et"ai] (|2002l . l2008al lbl) have de- 
veloped a formalism for handling radiative transfer in 
clumpy media and have applied it to AGN. Their model 
is able to reproduce some of the observed SED features, if 
a small number of clouds along equatorial rays is assumed. 
In particular, the infrared SEDs and strengths of the sili- 
cate emission features are now predicted to be largely inde- 
pendent of orientation and a large range of temperatures 
can coexist at a given distance from the central source, 
as re quired by new interferometric results (|Tristram et al.l 
|2007| . 2009). Nevertheless, the situation is far from being re- 
solved. Observations with the Spitzer Space Telescope have 



2 H. Landt, C. L. Buchanan and P. Barmby 



now started to show strong silicate emission features in 
quasars and absorption troughs in sou rces viewed edge-on, 
just a s predicted by Pier fc Krolik (e .g., SiebcnniorEcn ct al. 

2005; Wc cdman ot al 
,2006; Ogle et al.ll200e 



20051 : iHaas et all boOSl: iHao et al 



Clearv et all l2007l : iHao et al.| 



2005l:lBuchanan et al.ll2006l:IShi et al 



2007; 



Thompson et al.ll20 09: Hin er et al. 



Schweitzer et al.l 



2008 



2008 



Unravelling the structure of the dust obscuration is es- 
sential to our understanding of AGN. In this paper we use 
for the first time both photometric and spectroscopic Spitzer 
Space Telescope observations of radio-loud AGN, for which 
actual viewing angles can be determined, to investigate this 
topic. In Section 2 we introduce the sample and discuss the 
method used to derive viewing angles. The infrared data are 
presented in Section 3, whereas in Section 4 we confront the- 
ory with observations. Finally, in Section 5 we conclude and 
summarize our main results. 

Throughout this paper we have assumed cosmological 
parameters Ho = 70 km s~^ Mpc~^, Qm = 0.3, and f^A = 
0.7. 



2 THE SAMPLE 

The presence of strong radio jets in radio-loud AGN gives 
us a unique opportunity to derive actual viewing angles, 
8 (defined as the angle between the jet an d the observer's 
line o f sight), for this object class. Following iGhisellini et al.l 
1 19931 . we use here the method that combines apparent jet 
speeds, /3a (in units of the speed of light), with relativistic 
Doppler factors, S, to calculate viewing angles. 

If Pa is available from proper motion measurements and 
S can be separately estimated, then the viewing angle is: 



tan( 



2/3a 



13a' +S''-1 



(1) 



The largest uncertainty of this approach lies with the de- 
termination of 5, which involves several observables. Using 
the classical condition that the (predicted) synchrotron self- 
Compton (SSC) flux should not ex ceed the observ ed high- 
energy flux, 5 can be calculated as (|Marscheiil 19871 ): 



5 = f{a)Frr 



ln{ub/iym) 



F^& 



6+4a, „, 5-I-3Q 



l/(4+2a) 



(1 + ^) 



(2) 



where a is the spectral index of the thin synchrotron emis- 
sion (assumed to be 0.75), I'b is the synchrotron high- 
frequency cut-off (assumed to be 10^^ Hz), the function 
/(a) ~ 0.08q -I- 0.14 (assumed to be 0.2), and the red- 
shift z is known. The synchrotron flux F^ (in Jy) at the 
self- absorption frequency v^n (in Hz) of the core component 
with angular size Q^ (in mas) need to be determined from 
appropriate Very Large Baseline Interferometry (VLBI) ob- 
servations. Observables are also the high-frequency (X-ray) 
flux Fx (in Jy) at the (X-ray) frequency v^ (in keV) of this 
component. 

The most critical parameters in eq. ([2]) are the VLBI 
observables u^n and 0d, and these are in general difficult 
and expensive to determine. The determination of i/ni ide- 
ally requires a VLBI spectrum (i.e., high spatial resolution 
observations at different radio frequencies) of the strongest 
(core) component or a sophisticated spectral decomposition 
of a low-resolution spectrum to isolate the different (VLBI) 



components. On the other hand, observationally 0d will de- 
pend on the frequency itself (the higher the frequency, the 
higher the spatial resolution) and, therefore, will be an up- 
per limit only. This then makes 5 strictly speaking a lower 
limit a nd so the viewin g angl e an upper limit. 

As lGhiseUini et alj (|l993l ') (and also lRokaki et al.ll2003l ') 
have shown, this method samples a wide range of viewing 
angles, includ ing the regime of 20 ° < 6 < 60° that we are 
interested in. Vermeulen fc CohenI (| 19941 ) presented a large 
compilation of sources with measured jet proper motions. 
We have updated their list with later literature and have 
calculated viewing angles for all the sources. From the up- 
dated list, and restricting ourselves to 2; < 1 in order to keep 
the emitted far-infrared accessible, we have chosen for ob- 
servations with the Spitzer Space Telescope all sources with 
viewing angles 9 > 20° (12 objects; see Table [T|. We note 
that the source 3C 207 no longer obeys this criterion based 
on most recent proper motion data, but is kept in the sample 
for comparison. The cut in viewing angle avoids sources with 
infrared emission dominated by the relativistically beamed 
jet, so allowing for a reliable determination of the dust SED. 
All selected sources are radio quasars (i.e., type 1 AGN), ex- 
cept for 3C 84 {— NGC 1275), which is a radio galaxy (i.e., 
type 2 AGN). 

The total errors in our viewing angles are difficult to de- 
termine, since the errors in the Doppler factors 5 axe largely 
unconstrained. But, given the errors in the proper motions 
(Table [ij column (8)), we derive an average l ower limit on 
them of -^ 50%. In this respect, we note that iRokaki et al.l 
(2003) argue that it is unlikely that Doppler factors esti- 
mated using eq. ([2]) are wrong by a very large factor, since 
they observe excellent correlations between 5 and several 
emission line properties. We concur with this conclusion 
given our results in Section [4. 1.1 1 In addition, a consistency 
check that we carry out in Section 14.1.31 supports the no- 
tion that the trend in viewing angle that we obtain for our 
sample is qualitatively correct. 



3 OBSERVATIONS AND DATA REDUCTION 

We observed our sample with the Spitzer Space Telescope 
in Cycle 1 (ID: 3551). We used all its instruments and im- 
aged the sources with the Infrared Array Camera (IRAC; 
iFazio et al.ll2004l ) in four bands (3.6, 4.5, 5.8, and 8.0 ^m) 
and wit h the Mult i band Imaging Photometer for Spitzer 
(MIPS: [Rieke et aLllJOOJ ) in three bands (24, 70, and 160 
^m). In addition we obtained low-re solution spectrosco py 
with the Infrared Spectrograph (IRS; iHouck et al.ll2004l ) in 
staring mode using the appropriate short and long order 
modules in order to cover well the rest-frame wavelength 
region around 10 /im. 

In Tables [2] and O we list the details for the photomet- 
ric and spectroscopic observations, respectively. Some of the 
approved observations were allocated to other programs and 
were accessible to us only after the proprietary period ex- 
pired. We also note that we now include in this study all 
low-resolution (staring mode) IRS spectra available in the 
archive for our sources and not only those initially requested. 
In the following we describe the data reduction process and 
measurements. 



Table 1. General Properties of the Sample 



Object Name 


lAU Name 


R.A.(J2000) 


Decl.(J2000) 


z 


e 

[deg] 


[mas/yr] 


[mas/yr] 


Ref. 


Pa 


I'm 

[GHz] 


[Jyl 


ed 

[mas] 


Ref. 


[MJy] 


Ref. 


5 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 



III Zw 2* 
3C 47 
3C 84 

S5 0615+820 
3C 207 
3C 245 
3C 263 
3C 334 
3C 336 
4C +34.47 
4C +28.45 
3C 390.3 



0007+106 
0133+207 
0316+413 
0615+820 
0838+133 
1040+123 
1137+660 
1618+177 
1622+238 
1721+343 
1830+285 
1845+797 



00 10 31.0 

01 36 24.4 
03 19 48.1 
06 26 03.0 
08 40 47.6 

10 42 44.6 

11 39 57.0 
16 20 21.8 

16 24 39.1 

17 23 20.8 

18 32 50.2 
18 42 09.0 



+ 10 58 30 
+20 57 27 
+41 30 42 
+82 02 26 
+13 12 24 
+12 03 31 
+65 47 49 
+17 36 24 
+23 45 12 
+34 17 58 
+28 33 36 
+79 46 17 



0.090 
0.425 
0.018 
0.710 
0.681 
1.029 
0.646 
0.555 
0.927 
0.206 
0.594 
0.056 



28 
18 
31 
49 
8 
18 
48 
34 
23 
30 
24 
48 



0.166 
0.24 
0.266 
0.053 
0.334 
0.11 
0.06 
0.10 
0.10 
0.28 
0.13 
0.605 



0.016 
0.05 
0.050 
0.020 
0.030 
0.05 
0.02 
0.03 
0.07 
0.03 
0.04 
0.010 



MOJAVE 

VC94 

MOJAVE 

MOJAVE 

MOJAVE 

VC94 

VC94 

VC94 

H02 

VC94 

VC94 

MOJAVE 



0.99 
6.21 
0.32 
2.12 
12.94 
5.87 
2.23 
3.27 
4.94 
3.72 
4.50 
2.26 



15 
8 

28 
2 
8 
11 
8 
11 
8 
11 
5 
11 



0.79 
0.07 
6.50 
0.86 
0.60 
0.59 
0.12 
0.09 
0.01 
0.11 
0.50 
0.41 



0.15 
0.20 
0.35 
3.50 
0.28 
0.33 
0.37 
0.20 
0.20 
0.24 
0.30 
0.59 



B05 
H02 
R83 
FOO 
H02 
H87 
H02 
H92 
H02 
Ho92a 
Ho92b 
L81 



1.101 
0.363 
2.400 
0.040 
0.120 
0.160 
0.290 
0.342 
3.700 
1.467 
0.309 
7.600 



Pic05 
Ha06 
E06 
S97 
G03 
G03 
Ha02 
Ha99 
P96 
Pa04 
B97 
E06 



1.9 
0.4 
1.4 
0.4 
3.0 
1.7 
0.3 
0.4 
0.1 
0.2 
0.8 
0.1 



The columns are: (1) object name; (2) lAU name; (3) and (4) position, and (5) redshift from the NASA/IPAC Extragalactic Database (NED); (6) jet viewing angle calculated with 
the synchrotron self-Compton (SSC) formalism using the following observable s: (7) proper motion and (8) la error on p roper m otion, taken from reference in (9), whe re MOJAVE: 
large project for 'Monitoring Of Jets in AGN with VLBA Experiments' (e.e.. lKellermann et aLll2004l : iLister et aLlbOOQl) . VC94: IVermeulen fc CohenI {TgH), and H02: lHough et al.l 
1I2OO2I) : (10) apparent jet spee d calculated from colum n (7); (11) synchrotron self- absorp tion frequency; (12) f lux an d (13) size at (11) of the core c omponent; (1 4 ) refe rence for 
columns (11)-(13 ) , wher e B05: iBrunthaler et al.l ||2005^. F OO: iFey &: Chariot 1I2OOOI) . H87: 'Hough & Readhead! lll987f) . H92: iHough et al] l l 19921) . H02: iHough et all 1 I2OO2I) . Ho92a: 
Hooimever et al.l lll992d) . Ho9 2b: iHooimever et al.l lll992lj). L81:lLinfieldl lll98j) . R83: Readhead et al. (1983); (15) X-ray fl ux at 1 keV, taken from the refe rence in (16), where B97 : 
Brinkmann et al. Nl99'7l) . E06: lEvans et al.| ||200^ , G03: [Ga mbill et al.' (2003*), Ha99: Hardcastlc & WorralJ 119991) . Ha02: iHardcastle et akl liooi), Ha06: iHardcastle et ahl 1I2OO6I ). 

S97: .Sambruna 11997.) ; and (17) relativistic Doppler factor calculated from columns (11)-(13) and (15). 



P96: Pricto (199^), Pa04: |Page et al.l 1 I2OOI) . Pic05: |Picoiicem et al.. (.2005 ). 
position of the northern most source in a galaxy triplet 
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Table 2. Spitzer Photometry Journal of Observations 



Object Name 


ID 


IRAC 
observation 


texp. 




24 Aim 






MIPS 
70 /xm 






160 fim 








date 


[sec] 


ID 


observation 
date 


tcxp. 
[sec] 


ID 


observation 
date 


texp. 

[sec] 


ID 


observation 
date 


texp. 
[sec] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


III Zw 2 


3551 


2004 Dec 15 


4x 2 


86 


2004 Dec 26 


2x 3 


86 


2004 Dec 26 


2x 3 


3551 


2005 Aug 2 


10x10 


3C47 


3551 


2005 Jan 16 


4x 2 


3551 


2005 Jan 29 


3x 3 


3551 


2005 Jan 29 


3x10 


3551 


2005 Jan 29 


10x10 


3C 84 


3228 


2005 Feb 20 


5x30 


3551 


2005 Feb 26 


3x 3 


3551 


2005 Feb 26 


3x 3 


3418 


2005 Feb 1 


5x10 


S5 0615+820 


3551 


2004 Nov 20 


4x12 


3551 


2005 Mar 3 


3x 3 


3551 


2005 Mar 3 


3x10 


3551 


2005 Mar 3 


10x10 


3C 207 


3551 


2005 May 10 


4x12 


74 


2005 Apr 10 


1x10 


74 


2005 Apr 10 


4x10 


74 


2005 Apr 10 


4x10 


3C 245 


40072 


2008 Jun 10 


1x30 


40072 


2008 Jan 6 


1x10 














3C 263 


3551 


2004 Nov 1 


4x 2 


74 


2005 Apr 9 


Ix 3 


74 


2005 Apr 9 


1x10 


74 


2005 Apr 9 


4x 3 


3C 334 


3551 


2005 Mar 26 


4x 2 


74 


2005 Apr 7 


1x10 


74 


2005 Apr 7 


1x10 


74 


2005 Apr 7 


4x10 


3C 336 


3551 


2005 Mar 27 


4x12 


74 


2005 Apr 12 


1x10 


74 


2005 Apr 12 


6x10 


74 


2005 Apr 12 


4x10 


4C +34.47 


3551 


2005 Mar 30 


4x 2 


3551 


2005 Apr 6 


3x 3 


3551 


2005 Apr 6 


3x10 


3551 


2005 Apr 6 


10x10 


4C +28.45 


3551 


2004 Oct 8 


4x 2 


3551 


2004 Oct 17 


3x 3 


3551 


2004 Oct 17 


3x10 


3551 


2004 Oct 17 


10x10 


3C 390.3 


50763 


2008 May 13 


9x12 


3327 


2004 Oct 17 


4x10 


3551 


2004 Oct 18 


3x10 


3551 


2004 Oct 18 


10x10 



The columns are: (1) object name; for photometry with the Infrared Array Camera (IRAC) (2) program number, (3) observation date and (4) 
exposure time; for photometry with the Multiband Imaging Photometer for Spitzer (MIPS) (5) program number, (6) observation date, and (7) 
exposure time in the 24 fim band, (8) program number, (9) observation date, and (10) exposure time in the 70 ^m band, and (11) program number, 
(12) observation date, and (13) exposure time in the 160 fim band. 



3.1 The photometry 

The IRAC data were processed with the Spitzer Science Cen- 
ter (SSC) pipehne vl4.0, except for 3C 84 and 3C 390.3, for 
which we used the pipeline vl8.7.0 since it corrects arti- 
facts near bright sources. We measured flux densities from 
3.6 to 8.0 /im on the 'post-basic cahbrated data (BCD)' 
mosaics produced by the pipeline. For 3C 84 and 3C 390.3, 
mosaics made from the short-exposure images taken in high- 
dynamic range mode were used. We used aperture photom- 
etry and adopted a 10-pixel (12 arcsec) source radius, with 
background counts estimated in a 10-20 pixel (12-24 arcsec) 
radius annulus and subtracted. The sources 3C 336 and S5 
0615+820 had other objects nearby, therefore, we measured 
their fluxes in smaller 5-pixel (6 arcsec) apertures. 

We note that the 3.6 and 4.5 /im images of III Zw 2 show 
that this object may be a blend of two very close sources, 
and so flux densities at these wavelengths could be overes- 
timated. The 3.6 and 4.5 fim images of 3C 84 probably are 
contaminated by an extended component since no separate 
central point source is visible. For the source 3C 390.3 two 
data sets were available separated by about four years and 
we have considered both. However, the observations contem- 
poraneous with the MIPS photometry and IRS spectroscopy 
were severely affected by saturation, and, therefore, not use- 
ful for our purpose. 

All IRAC flux densities were aperture-corrected using 
the values in the IRAC Data Handbook (v3.0, 2006). No 
'array- location-dependent' photometric correction was ap- 
plied, since these are red sources and do not require such a 
correction. Table |4] lists our results. All sources have been 
detected in all four IRAC bands. We give la uncertainties 
estimated using the standard IRAF photometry formula. We 
note that these values do not include the err or in absolute 
calibration (a few per cent: lReach et al.ll2005l ). 

The MIPS data were reprocessed with the SSC pipeline 
vl6.1.0. Measurements on 24 fira images followed standard 
procedures using the post-BCD pipeline mosaics. We mea- 



sured flux densities in the standard aperture (35" radius) 
with the recommended aperture correction of 1.082 applied, 
except in the case of 3C 336 where a 7" radius aperture 
with a correction factor of 1.61 were used. The flux den- 
sity measured for 3C 84 (2.9 Jy) is formally above the sat- 
uration limit for 3-sec. exposures, but the object is only 
'soft-saturated' and the pipeline correctly replaced the sat- 
urated pixels with values from the 0.5-sec. exposures taken 
as part of the photometry AOR. Table 3] lists the results. 
All sources have been detected. We have computed uncer- 
tainties using the same method as for IRAC, which again 
do not include the error s in the absolute calibration (~ 2%; 
lEngelbracht et al.ll2007l ). 

For the measurements at 70 nm we used mostly the 
pipeline-produced 'filtered' mosaics, although in a few cases 
we remade the mosaics using time- and column-filtering 
of the BCDs to remove negative sidelobes around bright 
sources. Per the MIPS Data Handbook (v3.3.1) we mea- 
sured the flux densities in 35" radius apertures with sky 
annuli of 39-60 arcsec and an aperture co rrection of 1.22. 
Our r esults are listed in Table [l] Following ICarpenter et al.l 
((2009), we calculated uncertainties in the flux density using 
a similar equation to that for IRAC and MIPS-24, but ne- 
glecting Poisson noise and applying instead multiplicative 
factors to account for the noise correlation between pixels 
due to resampling during mosaicing and excess sky noise 
due to the data-taking proc edure. The absolute calibration 
uncertainty at 70 /im is 5% ({Gordon et al.ll2007l ). again not 
included in the tabulated values. The source 3C 84 is very 
bright at 70 /im but, as at 24 /im, is just below the satura- 
tion limit. Five sources were not detected at 70 /xm and we 
give only 3a upper limits. One source (3C 245) did not have 
MIPS- 70 data. 

For the measurements at 160 /xm we used the pipeline- 
produced 'flltered ' BCDs and made mosa ics using the 
MOPEX software (|Makovoz fc Marleauir2005l ) . However, for 
the bright source 3C 84 and the source 3C 390.3 that has 
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another bright object nearby we produced mosaics from the 
unfiltered images. Following the MIPS Data Handbook, we 
measured flux densities in 48" apertures with sky annuli of 
64-128 arcsec and an aperture correction of 1.60. The only 
exception was the source 3C 390.3. Due to the presence of 
a nearby brighter source separated by ~ 1 arcmin, we mea- 
sured its flux density instead in an aperture of 32" with an 
aperture correction of 1.97. Our results are listed in Table 
m Uncertainties for all sources were calculated as described 
above for the MIPS-70 data . The absolute calibrat ion un- 
certainty at 160 /iin is 12% (jStansberrv et al.ll2007^ and is 
not included in the tabulated values. The source 3C 84 is 
very bright at 160 /^m and at the saturation limit, therefore, 
its flux density is highl y uncertain and is taken to be 20% 
IStansberrv et al.ll2007l ). Seven sources were not detected at 
160 /^m and we give only 3(t upper limits. One source (3C 
245) did not have MIPS- 160 data. 

The archival MIPS data we used for four objects, 
namely, 3C 207, 3C 263 3C334 and 3C 336, were anal- 
ysed also by IClearv et al.l (|2007l l . Our flux densities in the 
24 pm , 70 /xm and 160 /im bands are consistent with theirs 
within 2(7, however, contrary to these authors we regard the 
sources 3C 207 and 3C 263 as undetected in the MIPS-70 
images and the sources 3C 334 and 3C 336 as detected in 
the MIPS-160 images. 



3.2 The spectroscopy 

The raw data were processed through the SSC pipeline vl5.3 
for the short-low (SL) modules and vl7.2 for the long-low 
(LL) modules. The BCDs produced by the pipeline were 
cleaned for rogue pixels using the IRSCLEAN softwarqj, 
then spectra were extracted using the SMART software 
IjHigdon et al.l [200J) . Where multiple exposures were ob- 
tained, the two-dimensional spectral images were median- 
combined. 

Two positions (nods) were observed for each target, and 
sky subtraction was performed on the images, by subtract- 
ing the off-source nod position image from the on-source im- 
age for each module. Spectra were extracted and the fluxes 
calibrated using the default tapered column apertures for 
each module. The modules at each nod position were merged 
and the edges and overlapping regions trimmed. The spec- 
tra from the two nod positions were then averaged to pro- 
duce the final spectrum. Uncertainty images provided by the 
pipeline and propagated through SMART were used to pro- 
duce the uncertainties on the final spectra. Our results are 
shown in Fig. [T] 



4 THE DUST GEOMETRY 

Fig. [T] shows the infrared SEDs for our sources plotted to 
have the same dynamic range. In the near-infrared we have 
supplemented the Spitzer data with ground-based photome- 
try in the J, H a nd Ks bands from the Two Micron All-Sky 
Survey (2MASS: ISkrutskie et"al]|2006l ). AU sources but one 



1 IRSCLEAN was written by the IRS GTO team (G. Sloan, D. 
Devost, &: B. Sargent). It is distributed by the Spitzer Science 
Center at Caltech. 



are included in the Point Source Catalogue of this survey 
and we list the flux densities in Table [5] For the undetected 
source S5 0615-1-820 we have derived 3cr upper limits based 
on the image plate specifications. 



4.1 General trends 

Before testing detailed theoretical models with our data 
(Section l4.2|l we want to first quantify the observations, with 
particular emphasis on revealing trends that any acceptable 
model will have to account for. 



4-1.1 The shape and width of the SED 

In the simplest approach we can approximate our obser- 
vations with a set of blackbodies. For this purpose we 
hav e fitted our data with the C routine MPFIT (version 
1.1: lMarkwardtll2009l ) . which uses the Levenberg-Marquardt 
technique to solve the least-squares problem. We have tried 
several numbers of blackbodies and have found that three 
or four components were required to obtain an acceptable 
fit. We have fitted for the temperatures and flux scalings of 
the individual blackbody components, meaning that our fits 
have either six or eight free parameters. Our results are listed 
in Table |S] and shown in Fig. [T] where dotted black curves 
indicate the individual components and the solid black curve 
their sum. 

We have included in the fit only the continuum part 
of the IRS spectrum, i.e., we have excluded strong narrow 
emission lines and the silicate emission features, and we have 
rebinned it to A log A = 0.05 /im in order to ensure a similar 
weighting between spectroscopy and photometry. However, 
in the case of the source 3C 245, for which neither MIPS-70 
nor MIPS-160 data were available, we have left the spectrum 
unbinned. We have treated all photometry upper limits as 
detections and have assumed their la values as the error. We 
have not considered the 2MASS photometric data points, 
since not only they appear to sample a different component 
(most likely the onset of the accretion disc) but they were 
obtained several years before our Spitzer observations. The 
two exceptions were III Zw 2 and 4C -1-34.47, for which we 
included the 2MASS Kb point, since it connected smoothly 
to the IRAC photometry and thus presented an important 
constraint on the hottest blackbody component. The IRAC 
data we use for the source 3C 390.3 was taken four years 
after the MIPS and IRS observations and variability by a 
factor of ~ 1.5 is observed. Therefore, assuming that vari- 
ability does not change the spectral slope, we have scaled 
these data to the IRS spectrum before including them in 
the fit. 

Our sample is evenly split into sources best-fit by three 
and four blackbodies. For the source 3C 245, the far-infrared 
part of the SED is not constrained and, therefore, we can- 
not exclude that four instead of three blackbodies might be 
required. In all sources but two the resulting fit is 'good' 
in a statistical sense (x^ of a few). In the two brightest 
sources (3C 84 and 3C 390.3) the resulting fit appears good 
to the eye and the large x?-values could be mainly due to the 
much smaller relative measurement errors involved. Result- 
ing temperatures for the hot, warm, cold and cool black- 
body components are in the ranges Thot '~ 1200 — 2000 
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Table 3. Spitzer Spectroscopy Journal of Observations 



Object Name 


ID 


observation 


peak-up 


SL2 


SLl 


ID 


observation 


peak-up 


LL2 


LLl 






date 




tcxp. 

[sec] 


tcxp. 

[sec] 




date 




texp. 
[sec] 


texp. 

[sec] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


III Zw 2 


86 


2005 Jul 10 


blue 


4x 30 


2x 60 












3C 47 


3551 


2005 Jan 11 


none 


- 


4x 14 


3551 


2005 Jan 11 


none 


4x120 


- 


3C 84 


14 


2004 Aug 30 


blue 


2x 14 


2x 14 


14 


2004 Aug 30 


blue 


4x 6 


4x 6 


S5 0615+820 


3551 


2004 Oct 23 


none 


- 


4x 60 


3551 


2004 Oct 23 


none 


4x120 


- 


3C 207 


3551 


2004 Apr 15 


none 


- 


4x 60 


74 


2005 Nov 15 


blue 


4x120 


4x120 


3C 245 












40314 


2008 Jan 10 


none 


- 


12x120 


3C 263 


3551 


2004 Nov 16 


none 


- 


4x 60 


74 


2004 Apr 14 


blue 


4x120 


4x120 


3C 334 


3551 


2005 Mar 14 


none 


- 


4x 14 


74 


2004 Jul 16 


blue 


4x120 


4x120 


3C 336 












74 


2004 Mar 22 


blue 


4x120 


4x120 


4C +34.47 


3551 


2005 Mar 14 


none 


- 


4x 6 


3551 


2005 Mar 14 


none 


4x 14 


- 


40 +28.45 


3551 


2004 Oct 23 


none 


- 


4x 60 


3551 


2004 Oct 23 


none 


4x120 


- 


30 390.3 


82 


2004 Mar 24 


blue 


1x240 


1x240 


82 


2004 Mar 24 


blue 


4x120 


4x120 



The columns are: (1) object name; for low-resolution spectroscopy with the Infrared Spectrograph (IRS) short order modules 
(2) program number, (3) observation date, (4) peak-up imaging array and exposure times for the slits (5) SL2 (5.2-8.7 ^m) and 
(6) SLl (7.4-14.5 ^m); for low-resolution spectroscopy with the IRS long order modules (7) program number, (8) observation 
date, (9) peak-up imaging array and exposure times for the slits (10) LL2 (14.0-21.3 fim), and (11) LLl (19.5-38.0 fim). 

Table 4. Spitzer Photometry Results 



Object Name 


F3.6 


0-3.6 


F4.5 


0-4.5 


F5.8 


05.8 


Fs.o 


os.o 


F24 


024 


F70 


070 


F160 


0-160 




[mjy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


III Zw 2 


22.2 


0.1 


28.1 


0.1 


35.3 


0.4 


48.8 


0.2 


138 


1 


128 


17 


59 


6 


3C 47 


5.3 


0.1 


6.9 


0.1 


9.2 


0.3 


11.2 


0.2 


36.8 


0.4 


37 


10 


<27 




30 84 


75.81 


0.09 


89.12 


0.06 


132.2 


0.2 


291.1 


0.2 


2930 


3 


3990 


33 


4350 


870 


S5 0615+820 


0.63 


0.01 


0.78 


0.01 


1.00 


0.05 


1.24 


0.02 


4.8 


0.3 


<18 




<21 




30 207 


1.54 


0.01 


2.09 


0.02 


2.75 


0.07 


3.94 


0.07 


13.2 


0.4 


<26 




<36 




30 245 


1.54 


0.01 


2.12 


0.01 


3.66 


0.03 


5.55 


0.03 


20.9 


0.6 










30 263 


5.5 


0.1 


7.3 


0.1 


8.9 


0.3 


11.5 


0.2 


26.7 


0.7 


<35 




<18 




30 334 


3.44 


0.09 


4.7 


0.1 


5.5 


0.3 


7.4 


0.2 


36.0 


0.3 


83 


11 


28 


7 


30 336 


0.76 


0.01 


1.04 


0.01 


1.44 


0.03 


1.93 


0.03 


3.9 


0.1 


19 


5 


62 


10 


40 +34.47 


9.0 


0.1 


11.9 


0.1 


13.8 


0.3 


18.2 


0.2 


59.9 


0.5 


<21 




<15 




40 +28.45 


2.19 


0.09 


3.4 


0.1 


3.8 


0.3 


5.4 


0.2 


18.2 


0.4 


<24 




<33 




30 390.3 


56.0 


0.1 


66.2 


0.1 


73.1 


0.2 


89.0 


0.2 


234.0 


0.2 


120 


6 


<24 





Table 5. Two Micron All-Sky Survey (2MASS) Fluxes 



Object Name 


J (1.235 fira) 


oj 


H (1.622 ^m) 


OH 


Ks (2.159 fim) 


"'Ks 




[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


III Zw 2 


4.5 


0.2 


6.3 


0.3 


12.9 


0.3 


30 47 


0.46 


0.05 


0.68 


0.07 


1.01 


0.09 


30 84 


12.4 


0.9 


16.6 


1.2 


20 


1 


S5 0615+820 


<0.13 




<0.22 




<0.26 




30 207 


0.34 


0.05 


0.42 


0.06 


0.64 


0.07 


30 245 


0.40 


0.05 


0.41 


0.08 


0.51 


0.08 


30 263 


1.87 


0.07 


1.79 


0.08 


2.30 


0.09 


30 334 


0.96 


0.05 


1.10 


0.08 


1.54 


0.08 


30 336 


0.40 


0.04 


0.33 


0.07 


0.45 


0.08 


40 +34.47 


2.72 


0.06 


3.03 


0.09 


4.7 


0.1 


40 +28.45 


0.32 


0.05 


0.53 


0.07 


0.73 


0.06 


30 390.3 


5.8 


0.2 


9.4 


0.3 


13.7 


0.4 
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Figure 1. Infrared spectral energy distributions plotted as rest- frame wavelength versus observed integrated flux. Filled black and green 
circles indicate Spitzer IRAC and MIPS and 2MASS photometry, respectively. The Spitzer IRS spectrum is shown in red, with its errors in 
yellow and the locations of the silicate features marked (red dotted lines). The dotted and solid black curves show the best-fit blackbodies 
and their sum. resnectivelv. The fit for the source .SC ."^QO-.S included the IR,AC data scaled to the IRS snectrum fonen black circles V 
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Table 6. Blackbody Fit Results 



Object Name* 


e 


That 


-'warm 


-'cool 


Tcold 


/hot 


/warm 


/cool 


/cold 


Xl/doi 




fdegl 


fKl 


[K] 


fKl 


fKl 


[erg/s/cm^] 


[erg/s/cm^] 


[erg/s/cm^l 


[erg/s/cm^] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


3C 207 


8 


1459 


509 


133 


- 


1.33e-12 


1.05e-12 


2.46e-12 


- 


5.3/8 


3C 47 


18 


1251 


456 


149 


- 


4.48e-12 


2.20C-12 


4.58e-12 


- 


1.62/7 


3C 245 


18 


1543 


501 


139 


- 


1.47C-12 


2.78C-12 


2.95e-12 


- 


8.1/98 


3C 336 


23 


2087 


842 


210 


58 


5.56e-13 


5.17e-13 


4.31e-13 


1.27e-12 


1.27/4 


4C +28.45 


24 


1175 


403 


146 


- 


2.27e-12 


8.92e-13 


2.78e-12 


- 


1.99/7 


III Zw 2 


28 


1561 


514 


176 


69 


1.87e-ll 


1.37e-ll 


1.37e-ll 


2.63e-12 


3.8/7 


4C +34.47 


30 


1702 


632 


189 


- 


6.80C-12 


5.27C-12 


6.53C-12 


- 


4.2/7 


3C 84 


31 


1290 


321 


138 


64 


6.65e-ll 


1.16e-10 


3.07e-10 


l.OOe-10 


1491/8 


30 334 


34 


1733 


799 


207 


98 


2.31e-12 


1.85e-12 


4.00e-12 


1.76C-12 


3.7/7 


30 263 


48 


1581 


666 


247 


133 


4.37e-12 


2.27C-12 


2.04e-12 


1.80e-12 


1.42/7 


30 390.3t 


48 


1274 


488 


183 


- 


3.04e-ll 


1.08e-ll 


2.53e-ll 


- 


757/11 


S5 0615+820 


49 


1654 


533 


179 


68 


5.22e-13 


3.13C-13 


5.83e-13 


5.36e-13 


0.55/6 



The columns are: (1) object name; (2) jet inclination angle; temperature of the (3) hot, (4) warm, (5) cool, and (6) cold 
blackbody component; peak flux of the (7) hot, (8) warm, (9) cool, and (10) cold blackbody component; and (11) reduced x^ 
value and number of degrees of freedom. 

* objects presented in order of increasing jet viewing angle 
1 fit performed to the IRAO data scaled to the IRS spectrum 




« 3C 207 
« III 7i* 2 



^-•••^»f«M(J?^'.^t 



> 3C 336 
• 4C +34.47 
« 3C 334 
« 3C 263 

> ,S5 0615-820 




1 10 100 

Figure 2. Infrared spectral energy distributions from Fig. [T]normalized at ~ 1 fim, grouped by the relative strength of the hot blackbody 
component and shown from left to right in decreasing order. The solid line indicates the sum of the best-fit blackbodies and has been 
omitted for the sources 30 207, 30 245, 40 +28.45, and S5 0615+820, for which the long-wavelength end is not well constrained. The 
locations of the silicate features are marked by the red dotted lines. 



K, Twarm ~ 300 - 800 K, Tcooi ~ 150 - 250 K, and 
Tcoid ~ 60 — 150 K, respectively. Note that the hottest 
blackbody component reaches values that are typical of the 
dust sublimation t emperature fo r most grain compositions 
(?a 1000 - 2000 K; ISalpeter|[l977l ). 

Three important trends are revealed by this simplistic 
approach. Firstly, the strengths of the individual blackbod- 
ies relative to each other vary substantially between sources, 
giving the impression that they indeed sample discrete com- 
ponents of a certain temperature rather than a single com- 
ponent with a smooth temperature distribution. Secondly, 
the cool (~ 200 K) blackbody component is prominent in 
all sources and appears roughly as strong or stronger than 
the hot blackbody component. Only in the source 3C 263 
is this behaviour reversed. And thirdly, the warm (~ 500 
K) blackbody component is often weak, thus introducing a 



sharp 'dip' in the SED just blueward of the 10 pm silicate 
feature. We also note that, although in most cases we did 
not include the 2MASS Ka point in the fit, it is often ap- 
proximated well by the hot blackbody component. 

In order to investigate if orientation determines to some 
degree the relative strengths of the individual blackbody 
components in a given source, we have normalized the SEDs 
at ~ 1 /xm and have grouped them in order of the relative 
strength of the hot blackbody component. In Fig. |2]we show 
our sources in three distinct groups with the average relative 
strength of the hot blackbody component decreasing from 
left to right. This representation of the SEDs suggests two 
important results. Firstly, no trend with orientation seems 
to be present for the hot dust. In particular, although the 
two sources with the strongest hot blackbody components 
(3C 47 and 3C 245) are among those with the smallest esti- 
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Figure 3. The ratio between the peak fluxes of the warm and 
hot blackbody components versus the jet viewing angle. 



relative to that of the hot component increase. This result is 
contrary to t he expectations of s mooth-density torus mod- 
els (see, e.g.. iPier fc Krolik|[l99^ . their Figs. 5 and 7). The 
source 3C 84 stands out in both Figs. [3] and |4] as having a 
relatively weak hot dust component. 

With the help of Fig. [5] we can identify differences and 
similarities between sources. Of special interest are possible 
shortcomings in the fits for those five sources that have two 
upper limits at the long-wavelength end, namely, 3C 207, 
4C -H28.45, 4C 4-34.47, 3C 263, and S5 0615-^820, and the 
poorly constrained source 3C 245. Firstly, we note that the 
values of the upper limits for the two sources 4C 4-34.47 
and 3C 263 are relatively low and, although located at dif- 
ferent rest-frame wavelengths, appear to sample a similar 
cool blackbody component. Therefore, they are unlikely to 
be far from the true values. Secondly, based on the similar- 
ity between the SEDs of the sources 3C 207, 4C 4-28.45 and 
S5 06154-820, and those of the sources III Zw 2, 3C 390.3, 
and 4C 4-34.47, respectively, the peak of the cool blackbody 
component of the former is unlikely to be overestimated by 
factors > 2. Finally, judging from a comparison between the 
SEDs of the sources 3C 245 and 3C 47, whereas the peaks of 
their cool blackbody components appear similar, the peak 
of the warm blackbody component of the former could be 
overestimated by a factor of '^ 2. 
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Figure 4. Same as Fig.[3]for the cool blackbody component. This 
component is likely to be overestimated for three sources (3C 207, 
4C -1-28.45, and S5 0615+820; encircled points). 



mated viewing angles, the main prediction of torus models 
assuming continuous density distributions, namely that the 
relative emission from h ot dust increases w ith decreasing 
viewing angle (see, e.g.. IPier fc Krolild Il992l . their Figs. 4 
and 5), is in general not observed for our sample. 

Secondly, it appears that at any viewing angle the SEDs 
can have widely different shapes. However, we find trends 
with orientation for the relative strengths of both the warm 
and cool dust components. These trends, which are illus- 
trated in Figs. [3] and [4l respectively, suggest that as the 
viewing angle decreases the warm and cool dust emissions 



4-1.2 The silicate features 

A characteristic property of interstellar dust are the spec- 
tral features centered at 9.7 /^m and 18 /iin due to silicates. 
In AGN, these features are expected to be produced by the 
dusty torus. A pertinent problem inherent to the measure- 
ment of silicate features remains the correct placement of the 
continuum, which ultimately decides if they are seen in emis- 
sion or in absorption. This task is still challenging because 
the available spectra rarely cover a wavelength region large 
enough to see a sizeable portion of the continuum around 
the features, which are relatively broad and often shifted 
in wavelength (see Fig. [T|. The advantage of our data set, 
however, is that the IRAC and MIPS photometry consider- 
ably extend the wavelength range of the IRS spectrum, thus 
allowing us to detect the overall continuum. In particular, 
our data shows that the 'dip' in the SED blueward of ~ 10 
^m observed in most sources, which could be interpreted as 
blueshifted silicate absorption, is in fact the result of a warm 
dust component that is relatively weak. 

Based on the overall continuum, we detect silicate emis- 
sion in all our sources. The only exception is the source 
3C 245, for which we do not detect strong silicate emission 
and interpret the IRS spectrum as being dominated by the 
warm and cool dust components. We have measured the 
properties of both the 10 /im and 18 pan silicate emission 
features after subtracting from the IRS spectrum the over- 
all continuum fitted in Section [4. 1.1 1 and removing superim- 
posed strong narrow emission lines. Our results are listed in 
Table [T] We detect the 18 /jm silicate feature in 4/8 sources, 
for which the IRS spectrum covers its location. In the re- 
maining four sources this feature appears to be swamped by 
the cool blackbody component, which, based on its temper- 
ature (~ 200 K), peaks around this wavelength. 

Our first noteworthy finding is that the 10 /xm sili- 
cate emission feature is considerably redshifted in all our 
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Table 7. Properties of the Silicate Emission Features 



Object Name* 


e 




Silicate 10 


fim 






Silicate 18 fj.m 








flux 


luminosity 


Wx 


center 


flux 


luminosity 


Wa 


center 






[erg/s/cm^] 


[erg/s] 


[Mm] 


[^m] 


[crg/s/cm^] 


[crg/s] 


[Mm] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


3C 207 


8 


4.42e-14 


8.91e+43 


0.65 


10.31 


- 


- 


- 


- 


3C 47 


18 


2.87e-13 


1.86e+44 


1.30 


10.65 


7 


7 


7 


7 


3C 336 


23 


8.45e-15 


3.63e+43 


0.37 


10.13 


- 


- 


- 


- 


4C +28.45 


24 


1.30e-13 


1.91e+44 


1.34 


10.61 


7 


7 


7 


7 


III Zw 2 


28 


2.38e-13 


4.79e+42 


0.16 


10.53 


7 


7 


7 


7 


4C +34.47 


30 


4.66e-13 


5.75e+43 


0.93 


10.44 


>1.04e-13 


>1.29e+43 


>0.29 


7 


3C 84 


31 


6.83e-12 


5.01e+42 


0.47 


10.83 


1.02e-ll 


7.41e+42 


0.61 


16.95 


3C 334 


34 


7.94e-14 


9.77e+43 


0.43 


10.53 


- 


- 


- 


- 


3C 263 


48 


1.05e-13 


1.86e+44 


0.56 


10.92 


4.31e-14 


7.76e+43 


0.34 


17.64 


3C 390.3t 


48 


1.22e-12 


9.12e+42 


0.59 


10.37 


4.13e-13 


3.09e+42 


0.27 


18.39 


S5 0615+820 


49 


>l.lle-14 


>2.51e+43 


>0.42 


>10.61 


7 


7 


7 


7 



The columns are: (1) object name; (2) jet viewing angle; for the 10 /im silicate emission feature (3) integrated flux, (4) 
luminosity, (5) rest-frame equivalent width, and (6) rest-frame center of the emission; for the 18 fj,m silicate emission 
feature (7) integrated flux, (8) luminosity, (9) rest-frame equivalent width and (10) rest-frame center of the emission. A 
question mark indicates that the feature is not covered by the spectrum. 

* objects presented in order of increasing jet viewing angle 

t continuum flt used the IRAC data scaled to the IRS spectrum 



sources. Instead of at the expected rest-frame wavelength 
of ~ 9.7 fira, we find its center to lie in the range of '^ 
10.1 - 10.9 fim (Table[71 column (6)). We note that strongly 
redshifted silicate emi ssion features were observed previ- 



ously in a few quasars (ISiebenmorgen et al.ll2005l : iHao et al.l 
I2OO5I : ISchweitzer et al.l I2OO8I ) , but never consistently in an 



entire sample. On the other hand, the center of the 18 /xm 
silicate emission feature in the three objects for which it can 
be determined is strongly shifted in only one source (3C 84) 
and in this case blueward. 

From circumstellar dust studies we know that the exact 
wavelength location of the silicate emission maximum de- 
pends on the grain size and composition, with larger grain 
sizes and larger amounts of crystalline dust (over amorphous 
dust) expect ed to enhance the em issivity at longer wave- 
lengths (e.g.. lBouwman et al.ll200in . In order to investigate 
if this assumption holds also for AGN, we have plotted in 
Fig. [5] the luminosity of the 10 ^m silicate feature versus 
its rest-frame center. Indeed, we find that stronger silicate 
emission is associated with a stronger redward displacement. 
However, the two lowest-luminosity sources (III Zw 2 and 
3C 84) do not appear to follow this trend. 

We note that in the current literature authors often 
study the (emission) equivalent widths of the silicate fea- 
tures rather than their luminosities. This approach needs 
to be taken with caution, since, by definition, the equiva- 
lent width depends strongly on the continuum fiux. And, 
as the SEDs in Fig. [T] show, the 10 pm silicate feature is 
located such that its continuum flux is given by the rela- 
tive strengths of the warm and cool blackbody components, 
which vary strongly between sources. In fact, the relation 
observed in Fig. [5] is not evident if we use rest-frame equiv- 
alent width values instead of luminosities, and our sample 
shows a trend for higher equivalent widths to be measured, 
the weaker the warm blackbody component is relative to the 
cool one. 
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Figure 5. The luminosity of the 10 fim silicate emission feature 
versus the feature center. Arrows indicate limits. 



4.1.3 The jet emission 

We have selected our sample of radio-loud AGN based on 
their relatively large jet viewing angles {8 > 20°). At these 
orientations the relativistic enhancement of the jet emis- 
sion is expected to be very low, ensuring that the observed 
infrared SED is dominated by emission from the putative 
dusty torus. In fact, based on the calculated relativistic 
Doppler factors S (Table[l] column (16)), we expect relativis- 
tic beaming of the integrated jet flux, which is proportional 
to (5'', by factors of < 10 in three sources (III Zw 2, 3C 84, 
and 3C 245) and relativistic debeaming for the remainder. 
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Figure 6. Infrared spectral energy distribution for 3C 207 from 
Fig.[l]extended to radio frequencie s. An extrapolation of the radio 
core spectr um [8 GHz VLBI data JHough et al.|[200 J ) and 5 GHz 
VLA data JHardcastle et aUl2004h : blue filled circles] to higher 
frequencies using a typical blazar SED (blue solid curve) shows 
that even in our strongest relativistically beamed source the jet 
does not dominate the infrared emission. 
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Figure 7. The ratio between the integrated core radio luminosity 
and the peak luminosity of the hot blackbody component versus 
the jet viewing angle. 



The exception is the source 3C 207, for which, based on 
improved proper motion data, we now calculate a relatively 
small inclination angle (6 ~ 8°). Its relativistic Doppler fac- 
tor is (5 ^ 3, meaning that the relativistic jet flux enhance- 
ment is expected to be a factor of ~ 80. However, as Fig. [6] 
shows, even in the strongest relativistically beamed source 



in our sample the infrared SED is dominated by thermal 
emission. An extrapolation of the core radio spectrum us- 
ing a typical bl azar SED that generally peaks at a few ^m 
iGiommi et al.iii2002 ') predicts mid-infrared fluxes a factor of 
~ 40 lower than observed. 

Further evidence that the infrared SEDs of all our 
sources are dominated by thermal rather than non-thermal 
emission comes from their similar, relatively bulgy appear- 
ance and the lack of strong, short-term variability. The IRS 
spectroscopy overlaps in wavelength with either the IRAC 
and/or the MIPS photometry, and an impressive consistency 
is evident between the two data sets, which were taken sev- 
eral months apart (see Tables [2] and |3]) . Variability is de- 
tected in the source 3C 390.3, however, not of the strength 
typical of blazars. The flux increase is a factor of < 2 over a 
period of four years. 

Given that the infrared emission in our sources is ther- 
mal and that the hot dust emission seems to be emitted 
isotropically (see Section|4]TTT]l, we can now carry out a con- 
sistency check for the calculate d viewing angles. Al o ng th e 
lines of argument presented by IWills fc BrothertorJ ([1995|) , 
we plot in Fig.[7]the ratio between the integrated core radio 
luminosity (calculated from the data listed in Table [ij and 
the peak luminosity of the hot blackbody component (cal- 
culated from the fluxes listed in Table (6] column (7)) versus 
the jet viewing angle. In such a diagram, we expect that the 
stronger a source is relativistically beamed, the higher its 
ratio between beamed and isotropic emission. Fig. [7] shows 
that the resulting trend in viewing angle for our sample 
is qualitatively correct; the smaller the viewing angle, the 
higher the ratio between (beamed) core radio power and 
(isotropic) dust luminosity. The only pronounced exception 
is the source S5 0615-1-820, for which the viewing angle ap- 
pears to be highly overestimated. 

4.2 Testing CLUMPY torus models 

Due to the shortcomings in AGN torus models invoking con- 
tinuous dust distributions, new models ba sed on clumpy me- 



dia h a ve rec e ntly been put forw a rd (e.g., Schartmann et al.l 



20051, l2008l: iHonig et all l2006l : iHonig fc Kishimotd l201o[ 



Nenkova et al.ll2002l . l2008al lbl). The most detailed of these 
models are those of Nenkova and collaborators (named 
CLUMPY), and they can be accessed on-linqj. In short, 
these authors solved the radiative transfer problem in 
clumpy media by assuming that the medium is composed 
of clouds that are individually optically thick (rv3>l), that 
each cloud can be considered a point source of intensity 5c, a, 
and that the cloud distribution obeys Poisson statistics. In 
this case, the escape probability of the emitted radiation can 
be approximated as Pose — e~^'^, with A^'t the total number 
of clouds along the line of sight, and the intensity at a given 
location becomes: 



I?(s) = 



-iVT{a,a') 



S,,xNc{s')ds', 



(3) 



where A'^c is the radial cloud density (i.e., the number of 
clouds per unit length). With this formalism the only dif- 
ference between the clumpy and continuous cases is that 

^ See jhttps : //newton ■ pa . uky . edu/ - clumpyweb/ 1 
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Figure 8. The 10 /im silicate emission feature and the spec- 
tral slope blueward of it are well approximated by redshifted 
CLUMPY models with standard ISM dust composition (code 
ISM-D; top panels, solid lines). Dust composed of silicates with 
various materials added (code OHM; bottom panels) does not 
reproduce the observations as well, independent of whether we 
apply a redshift (solid lines) or not (dashed lines). 



optical depth rv is replaced by its effective equivalent 
A't(1 — e~^^) ~ A'^T (for ry^l) and the absorption coef- 
ficient is replaced by A^c- The main challenge lies with the 
calculation of the clump source function, 5c, a, that needs 
to properly take into account the illumination profile of the 
individual clouds and the effects of cloud shadowing. 

In this section we wish to test the CLUMPY models, 
in particular how well they approximate the entire infrared 
SED, how unique a fitting set of parameters is and what 
physics can be extracted. For simplicity, we will limit our- 
selves to models invoking a torus geometry with a Gaussian 
angular cloud distribution, i.e., we will not consider a spher- 
ical distribution or the sharp-edge geometry. 



4-2.1 Constraining the parameters 

The CLUMPY torus models are available for two differ- 
ent dust compositions and have six parameters that can 
be adjusted: the optical depth of the cloud, ry, the aver- 
age number of clouds along an equatorial line of sight, A'^o, 
the power-law index of the radial cloud distribution, q, the 
angular width of the Gaussian cloud distribution, a, the ra- 
tio between the outer and inner torus radius, Y, and the 
viewing angle, 9. The dust composition is assumed to be a 
mix of silicates and graphite with silicates either as in the 
standard interstellar medium (IDraind 20031, code ISM- D) or 
with various materials added ( Ossenkopf et al.l Il992l . code 
OHM). 

In general, the problem that one faces when choosing a 
CLUMPY model that best-fits the data is that its parame- 
ters are highly degenerat e, with several permutations yield- 
ing very similar results (JAsensio Ramos fc Ramos Almeidal 




Figure 9. The height of the 10 fj,m silicate emission peak relative 
to the 'dip' blueward of it constrains best the number of clouds 
(top panels, where dotted, solid and dashed lines indicate A^o = 4, 
5 and 6, respectively), the power-law index of the radial cloud 
distribution (middle panels, where dotted, solid and dashed lines 
indicate g = 1, 2 and 3, respectively), and the angular width of the 
Gaussian cloud distribution (bottom panels, where dotted, solid 
and dashed lines indicate a = 30°, 45° and 60°, respectively). 



I2OO9I : iRamos Almeida et all I2OO9I : iNikutta et aD 120091 ) . In 
the following, we introduce an effective method to select a 
suitable model (by eye) and find that the high degeneracy is 
due to the fact that, with the exception of the dust compo- 
sition, all parameters constrain only the radial mass density 
profile (and so the radial temperature profile). This then 
means that, although CLUMPY cannot uniquely determine 
the AGN torus structure, it is a versatile tool with possible 
application to other astrophysical objects. 

As a first step, we find that our observations can con- 
strain the dust composition. As Fig.[8]shows, ISM-D models 
approximate well both the shape of the 10 pm silicate emis- 
sion feature and the spectral slope blueward of it, however, 
only if a redshift is applied (top panels). On the other hand, 
OHM models predict a narrower silicate feature than is ob- 
served and do not match well the blueward spectral fiux, 
independent of whether we apply a redshift or not (bottom 
panels). We have then used in the following only ISM-D 
models. In this respect, we note that our result diffe rs from, 
but does not contradict, that of lSirockv et al.l l|2008r ). These 
authors constrain the dust composition of ultraluminous in- 
frared galaxies based on the strengths of both the 10 ^m 
and 18 pm silicate features to be that of the OHM mod- 
els. However, a careful inspection of their Figs. 7-9 shows 
that their method cannot constrain the dust composition in 
type-1 AGN, such as our sources. 

Next, we find that the most convenient and reliable way 
to constrain the five geometrical CLUMPY parameters is to 
use the height of the 10 nm silicate emission peak relative 
to the 'dip' blueward of it. Then, this fiux ratio is larger, 
the smaller A^o, <1, "" and 6, and the larger Y are. The most 
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Table 8. CLUMPY Torus Model Fit Results 



Object Name* 


e 

[deg] 


^C 


TV 


No 


9 


-f'bol 
[erg/s] 


[pc] 


Mc 
[Mq] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 



3C 207 


8 


0.225 


10 


7 


3 


1.12e+46 


1.340 


8.5O+04 


3C 47 


18 


0.186 


10 


5 


2 


1.15e+46 


1.355 


1.3e+05 


3C 245 


18 


0.225 


150 


20 


3 


3.47e+46 


2.355 


l.le+07 


3C 336 


23 


0.202 


10 


5 


3 


9.77e+45 


1.250 


5.3e+04 


4C +28.45 


24 


0.167 


10 


4 


2 


1.41e+46 


1.503 


1.3e+05 


III Zw 2 


28 


0.230 


10 


7 


3 


1.48e+45 


0.486 


l.le+04 


4C +34.47 


30 


0.186 


10 


5 


2 


4.07e+45 


0.807 


4.7e+04 


3C 84 


31 


0.109 


150 


20 


3 


3.09e+44 


0.222 


l.Oe+05 


3C 334 


34 


0.067 


150 


20 


3 


9.33e+45 


1.222 


3.0O+06 


3C 263 


48 


0.230 


10 


5 


2 


3.55e+46 


2.383 


4.1e+05 


3C 390.3 


48 


0.175 


10 


5 


1 


9.55e+44 


0.391 


2.4e+04 


S5 0615+820 


49 


0.202 


10 


4 


1 


6.61e+45 


1.028 


1.3e+05 



The columns are: (1) object name; (2) jet viewing angle; for the CLUMPY component: 
(3) redshift, (4) cloud optical depth, (5) number of clouds along equatorial rays, (6) 
power-law index of the radial cloud distribution, (7) bolometric luminosity, (8) inner 
radius, and (9) mass in clouds (in solar masses). In all cases we assumed parameters 
0- = 45 and y = 10. 

* objects presented in order of increasing jet viewing angle 



pronounced differences are produced by changing A'^o, q, and 
a (Fig.O, and the strongest degeneracy seems to be present 
between A^o and 6 (i.e., we obtain a similar result, if we either 
increase A^o and decrease 9 or vice versa). Since all these pa- 
rameters effectively determine th e cloud number density (see 
eq. (2) of iNenkova et al.l l|2008bl )). and the aforementioned 
degeneracy serves to keep it constant, this then means that 
our flux ratio is larger, the smaller A^c- Once A^c is fixed, the 
selection of tv is straightforward, since models with small 
values predict much stronger 10 /xm silicate emission and 
much less radiation at far -IR wavelengths than models with 
large values (see Fig. 5 of iNenkova et al.ll2008bl ). 



4-2.2 Revealing the strengths and deficits 

Applying the CLUMPY models to the infrared SEDs of our 
sources, we find in all cases that they can approximate well 
the mid-IR part, but significantly underpredict the flux at 
both near-IR and far-IR wavelengths. Furthermore, it is al- 
ways necessary to apply a redshift to the models. We show 
the models (including the predicted AGN contribution) in 
Fig. [To] (black solid curves) and list the fltting set of param- 
eters in Table [S] 

The cold dust component unaccounted for by CLUMPY 
models is ubiquitously observed in AGN and has been as- 
sociated with dust from either starburs ts and/or the nar- 
row emission line regi on (NELR) (e.g. , iNetzer et al.l 120071 : 
ISchw citzcr et al. 2008; iMor et all [20091 ). The flrst alterna- 
tive is unlikely to hold for our sample for three reasons: 
(i) radio-loud AGN reside predominantly in luminous el- 
hpticals, which usually do not have signiflcant starbursts; 
(ii) we do not observe strong polycyclic aromatic hydrocar- 
bon (PAH) features in our IRS spectra, which are usually 
associated with starbursts; and (iii) the observed peak lu- 
minosities of the cold dust components are of the order of 
icoid ~ 10^" — 10^^ solar luminosities, implying that the star- 
burts would have to be unusually luminous. On the other 



hand, models for the dust emission from the NELR usually 
predict higher temperatures and, therefore, aim to account 
mai nly for the wavel ength region covered by CLUMPY mod- 
els (|Mor et al.ll2009h . 

Besides the cold dust component, CLUMPY cannot ac- 
count for the hottest dust. The reason for this might be sim- 
ply that its chemical composition differs from that assumed 
by the models. For example, it could be composed of pure 
graphite, which has a higher sublimation temperature than 
silicates and emits a spectrum close to a pure blackbody. 
Other authors have previousl y hinted at such an origin for 
the hot dust component (e.g.. lMor et al.ll2009l ). 

The redshifts required to displace the models are rel- 
atively large {zc ~ 0.07 — 0.23). Fig. 1111 shows these red- 
shifts compared to those we derive from the measured cen- 
ter of the 10 ^m silicate emission feature (Table O, with 
the latter redshifts calculated assuming the peak rest-frame 
wavelength of the ISM-D silicates (9.5 /im instead of 9.7 
/im). With the exception of the two sources that are mod- 
eled with a very high optical depth (3C 84 and 3C 334, see 
below), which makes the determination of zc rather uncer- 
tain, the CLUMPY redshifts are always larger. In order to 
understand this flnding we recall that we determine zc by 
adjusting the CLUMPY model to the location of the 'dip' 
blueward of the silicate feature and thus to the feature's 
blue wing, whereas the measurement of the emission cen- 
ter depends strongly on the width of the feature's top part. 
However, as CLUMPY shows (Fig. \^, the width of the 
top part is reduced as the cloud number density and so the 
absorption increases, which leads to an apparent blueshift 
of the emission peak and thus an underestimating of the 
(true) redshift. This radiative transfer effect then also offers 
an alternative explanation for the relation found in Fig. [5] 
Instead of a difference in grain size, it could be that the 
larger the observed redshift of the center, the less absorbed 
the silicate feature is, and so the larger its luminosity (for 
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Figure 10. Infrared spectral energy distributions from Fig. [T] with redshifted CLUMPY torus models overlaid (solid black curves). 
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the same AGN bolometric luminosity, as is the case for our 
sample; see be low) . 

Recently, iNikutta at al.1 ([2003) proposed that the ob- 
served redshifts of the 10 fim silicate emission features in 
AGN can be explained by an interplay between the radia- 
tive transfer effects illustrated in Fig. [12] and a rising con- 
tinuum underneath the feature. In particular, they showed 
that the appearance of the feature's top part can change 
from less to more peaked depending on if the region be- 
tween the 10 /im and 18 /im silicate features was included 
in the continuum fit or not. In this respect we note that we 
have included the middle continuum part in our fit and have 
measured the center (not the peak) of the feature, which is 
less sensitive to the continuum placement. Furthermore, as 
Fig. 1121 shows , the radiative transfer effects mentioned by 
iNikutta et al.l ((20091) decrease the observed redshifts and so 
cannot explain them (see also iHo nig ct al. 2010). 

The large majority of our sources are best modelled 
assuming a cloud optical depth of rv = 10, with only three 
objects (30 245, 3C 84, and 30 334) requiring a much larger 
value of TV = 150 in order to produce enough emission at 
large (A > 20 /xm) wavelengths. The number of clouds along 
equatorial rays and the index of the radial cloud distribution 
differ the most among our sources, with required values of 
A'o = 4,5,7 and 20, and g = 1,2 and 3, respectively. On 
the other hand, we find that it is not necessary to vary the 
ratio between the outer and inner torus radius and the width 
of the Gaussian distribution, with values of Y = 10 and 
a = 45° , respectively, giving satisfactory results. 

Several physical parameters can be derived from an ap- 
proximation of the infrared SED with OLUMPY models. 
The most important of these are the bolometric luminosity 
of the source, the inner radius of the OLUMPY component, 
and the observed mass in clouds (Table |S}. The bolometric 
luminosity can be readily obtained from the scaling of the 
model to the data, taking into account that the model is 
redshifted and thus the observed flux reduced by a factor of 
{1 + zc)^. The resulting values for our sources span only ~ 2 
orders of magnitude and we note that in particular the nine 
objects that follow a relation in Fig. [S] have similar values 
((logLboi) = 45.9 ±0. 1 erg s~^). Using the bolometric lu- 
minosity and eq. (1) of lNenkova et al.l (|2008a ) we have then 
calculated the inner radius of the CLUMPY component, Rin, 
assuming a dust sublimation temperature of Tsub — 1500 K. 
We obtain values in the range of -Rin ~ 0.2 — 2 pc. 

The calculation of the total mass in clouds follows eq. 
(6) of lNenkova et al.l (l2008ar ) and re lies on the cloud number 
density, A^c ■ iNenkova et al.l (|2008a ) give an analytic expres- 
sion for the mass in the case of the sharp-edge geometry that 
differs only slightly from the proper numerical integration in 
the case of a Gaussian toroidal geometry: 



Mc = muNu / Ac dV 



TV 2 

mn — -RinAo/q47r sin a, 

(TV 



(4) 



where ttih is the proton mass. Ah is the hydrogen col - 
umn density of a cloud, av=4.89e-22 cm^ JDrainel [20031 ) 
is the dust absorption cross-section in the V band, and 
7q = {Y^ - l)/(21nF), Y, and (21ny)/(l - y^) for the 
case of g = 1, 2, and 3, respectively. Since the degeneracy 
inherent to OLUMPY serves to keep the mass density con- 
stant, the determination of Mc is expected to be robust, i.e.. 
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Figure 11. The redshift of the center of the 10 fim silicate emis- 
sion feature versus the redshift required to displace the CLUMPY 
models. The solid line marks the locus of equality. Arrows indicate 
limits. 
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Figure 12. The emission peak of the 10 ^m silicate feature ap- 
pears increasingly blueshifted as self-absorption becomes impor- 
tant. In the examples shown the feature center has shifted from 
9.5 fira (top left panel) to 8.1 ^m (bottom right panel) (vertical 
dotted lines). 



one should get the same result with any fitting set of model 
parameters. 
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5 SUMMARY AND CONCLUSIONS 

We have presented infrared observations obtained with all 
three instruments on-board the Spitzer Space Telescope for 
12 radio-loud AGN, for which actual viewing angles can be 
determined. The results of our analysis of the infrared SED 
and the strength of the 10 /im silicate feature can be sum- 
marized as follows: 

1. The infrared SED is best-fit with a combination of 
three or four blackbodies. The resulting temperatures for 
the hot, warm, cool and cold components are in the ranges 
Thot ~ 1200-2000 K, Twarm ~ 300-800 K, T^ooi ~ 150-250 
K, and Tcoid ~ 60 — 150 K, respectively. 

2. We find trends between the emissions of the warm 
and cool dust components and the jet viewing angle indi- 
cating that the more the source is viewed face-on, the larger 
the amount of dust that comes into our line of sight. No 
such trend is present for the hot dust component. Both these 
results are contrary to the expectations of smooth-density 
torus models. 

3. Based on the overall continuum, we detect the 10 /xm 
silicate feature in emission in all our sources. The feature 
center is always observed to be strongly redshifted. 

4. We test the CLUMPY torus models of Nenkova and 
collaborators and find that they approximate well the mid- 
infrared part of the SED, but significantly underpredict the 
fiuxes at both near-IR and far-IR wavelengths. Furthermore, 
we find that the models can constrain the dust composition 
(in our case to that of the standard ISM), that they require 
relatively large redward displacements (^ 10% — 20% the 
speed of light) to match the observations, and that they are 
insensitive to the assumed geometry but give robust total 
mass estimates. 
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